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Q, The global potential energy surfaces of the first six singlet electronic states of 



< 



CO2, 1 — 3^yl' and 1 — 3^^" are constructed using high level ab initio calculations. 

In linear molecule, they correspond to X^T,^, l^A^, l^S^, and l^n^. The cal- 

' r-| ' 
O, culations accurately reproduce the known benchmarks for all states and establish 

C missing benchmarks for future calculations. The calculated states strongly interact 

CD 

r~| at avoided crossings and true intersections, both conical and glancing. Near degen- 

q 

ty^ eracies can be found for each pair of six states and many intersections involve more 

O 
*^ than two states. In particular, a fivefold intersection dominates the Franck-Condon 

>-. 

^S^ zone for the ultraviolet excitation from the ground electronic state. The seam of 



Ph 



this intersection traces out a closed loop. All states are diabatized, and a diabatic 



►^ 5x5 potential matrix is constructed, which can be used in quantum mechanical 

\^ calculations of the absorption spectrum of the five excited singlet valence states. 

O I- INTRODUCTION 

m 

^ This and the subsequentP paper (termed 'paper 11') describe the results of an ab initio 

^ quantum dynamical study of the absorption spectrum and the non-adiabatic dissociation 

^ mechanisms of carbon dioxide photoexcited with the ultraviolet (UV) light between 120 nm 

and 160 nm. A brief account of this work has already been published.'^' Paper II gives the 
motivation behind the study and discusses its main result — the quantum mechanical ab- 
sorption spectrum and its interpretation in terms of wave functions of metastable resonance 
states. The present paper sets the stage for paper II and describes the ab initio calculations 
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and the topography of the potential energy surfaces (PESs) involved in photodissociation. 
This information, which is often stashed away in supplementary online sections,'^ is a central 
and, indeed, an indispensable ingredient of a reliable dynamics calculation. The construction 
of ab initio PESs and their diabatization — which, without much exaggeration, amounts to 
learning the topography of PESs and their intersections by heart — is often more challenging 
than the subsequent quantum dynamical calculation. 

Photoabsorption from the ground electronic state X^H'^ of linear CO2 at wavelengths 
120 nm — 160 nm is due to the first five excited singlet valence states l^S^, l^H^, and 
1^ A^.^"- In the Cs group notation, appropriate for bent molecule, these states are 2, 3M' 
and 1, 2, 3^A". UV light excites CO2 into the region of multiple electronic degeneracies, nu- 
clear motion through which induces strong non-adiabatic couplings between electronic states. 
These couplings directly affect the observed absorption spectrum of the valence states and 
control distributions of photofragments over the final states. Their indirect influence appar- 
ently extends to shorter wavelengths where Rydberg transitions dominate: The combined 
experimental and theoretical analysis indicates^ that the manifold of coupled valence states 
acts as a 'sink' for the optically bright Rydberg states and affects their dissociation lifetimes. 

Electronic degeneracies in the Franck-Condon (FC) region, which have been the focus 
of several studies in the past,'2EEHlol g^j^g Qf ^^q types. Glancing intersections occur in the 
orbitally doubly degenerate l^Hg and l^A^ states which upon bending split into A' and A" 
components. Conical intersections (CIs) arise at the accidental l^Ug/l^Au crossing inside 
A' and A" symmetry blocks.'^'^EIISI j^ fact, CIs between valence states are ubiquitous and 
found far outside the near-linear FC region, at strongly bent geometries J ^ * ^ * ^° ' Together with 
local minima and saddles, these crossings are the principal features shaping the topography 
of the singlet valence states. 

The outline of the paper is the following. Section [IT] sketches the technical details of ab 
initio calculations. Next, the constructed adiabatic PESs are presented for the ground (Sect. 



Ill A) and the excited (Sect. IIIB) electronic states. Their intersections are discussed for 



near linear (Sect. IIIC) and bent geometries (Sect. |IIID), and put into perspective by a 



review of a network of closely spaced valence and Rydberg states (Sect. |IIIE). 



If the solution of the Schrodinger equation for nuclei is out of reach, the description of 



adiabatic surfaces in Sect. Ill would be the last step in theoretical ab initio analysis. If, on 



the other hand, one intends to treat nuclear dynamics quantum mechanically, ab initio PESs 



featuring CIs have to be diabatized.'^ This is especially desirable if — as in paper II — a 
discrete grid is used to represent the nuclear Hamiltonian, because the diabatized potential 
matrix is free from either divergent off-diagonal couplings or non-differentiable potential 
cusps. While general schemes for constructing approximate adiabatic-to-diabatic transfor- 
mations are established (see, for instance, Ref. [T2ll . their application to the valence states 
of CO2 is complicated by the number of CIs to be simultaneously treated. Simplifications 



are called for, as described in Sect. IV, in which the diabatic representation is constructed 



separately for bent CIs (Sect. IV A) and linear CIs (Sect. IV B). Section M concludes. 



II. ELECTRONIC STRUCTURE CALCULATIONS 

All ab initio calculations are carried out with the MOLPRO package.!^ The Gaussian 
atomic basis sets used in this work are due to Dunning.^ Previous studies indicate that 
diffuse functions should be added to the basis sets on oxygen and carbon atoms in order 
to account for the mixed valence-Rydberg character of the 11 state.'^'^*^ A series of tests 
was conducted, in which s, p, d etc. basis functions of triple and quadrupole zeta quality 
were selectively augmented by one or two diffuse functions. The pre-computed doubly 
augmented correlation consistent polarized valence quadrupole zeta (d-aug-cc-pVQZ) basis 
set,'^as implemented in MOLPRO, was found computationally most stable and selected for 
calculations of global PESs. 

Three-dimensional (3D) PESs of states 1,2,3M' and 1,2,3M" are calculated at the 
internally-contracted multireference configuration interaction singles and doubles (MRD- 
CI) level, based on state-averaged full-valence complete active space self-consistent field 
(CASSCF) calculations with 16 electrons in 12 active orbitals and 6 electrons in three fully 
optimized closed-shell inner orbitals. The electronic configuration of the ground state X^'E'^ 
is ([core]2cr^3cr^4cr|3cr^l7r^l7r^). The dominant electronic excitations, leading to the lowest 
excited states, include Ivr^ — )■ 2p2nl (giving states l^'^A^ and l^T,~) and Ivr^ — )■ SsScr^ (giving 
state l^IIg). Active orbitals in CASSCF comprised 2cru — 4cr„, Sag — ^Og.XiXu — 2vr„, and Ivr^. 
In Cg symmetry, used in the calculations, these are 4a' — 12a' and la" — 3a". In the MRD-CI 
step all 16 valence electrons were correlated. The maximum numbers of open shells allowed 
in the MRD-CI calculations were 8 in the reference space and 12 in the internal space. This 
lead to 38159928 contracted configurations. The Davidson correction was applied in order 



to account for higher-level excitations and size-extensivity.'^ 

Adiabatic energies are calculated on a 3D grid of the two C-0 bond lengths Ri^2 and 
the OCO bond angle aoco- R2 ^ [1.9ao, 3.0ao] (step size equals 0.1 oq), -Ri G [-R25 6.2ao] 
(step size varies between 0.1 oq and 0.4 oq), aoco ^ [70°, 179.9°] (step size varies between 
2° and 10°). Additionally, many cuts in the (i?i,i?2) plane are computed for angles aoco 
between 60° and 0° within the continuing effort to construct a balanced description of both 
CO + O and C + O2 arrangement channels. At present, the grid comprises 4800 symmetry 
distinguishable points. The resulting energies were scanned in one and two dimensions for 
obvious errors. The list of corrected adiabatic energies was subsequently interpolated using 
3D cubic splines and also used for constructing the quasi-diabatic representation. Missing 
energies for aoco < 70° in the dissociation channels were obtained from those for aoco > 70° 
using trigonometric extrapolation. 

Absolute intensity calculations of paper II require transition dipole moments (TDMs) 
with the ground state X. Components [iJ'x^fJ'y, fJ'z) of the TDM vector are calculated, for 
each electronic state, on a 3D grid Ri^2 = [1.9ao, 2.4ao] and a = [165° — 179°] covering 
the spot over which the vibrational ground state in X is delocalized. The molecular axes 
in these calculations are chosen such that x' is orthogonal to the molecular plane, z' runs 
along one of the CO bonds and y' _L z'. For A' states, the in-plane components are generally 
non-zero, while for A" states, it is the fix' component which carries the transition. 

III. PROPERTIES OF THE VALENCE PESs AND THEIR CROSSINGS 

A. Ground electronic state 

The ~ 7.5 eV deep adiabatic ground state PES supports three structural isomers: The 
familiar linear OCO molecule is the global equilibrium, while the carbene-like bent OCO 
and the linear COO are the two local ones. Table [T] summarizes the characteristic features 
of the X^Ti'f state at the three equilibria and compares them with the previous ab initio 
studies and with the available experimental data. 

The vicinity of the global minimum is of capital importance for the environmental 
chemistry.'^ The calculated equilibrium CO bond distance in linear OCO, Re = 2.1991 ag, 
agrees well with the experimental value of 2.1960 oq. The accuracy of the vibrational zero- 



point energy (ZPE) and the vibrational transitions frequencies is assessed in Table III] which 
compares energies of the low lying vibrationally excited states in rotating CO2 (the to- 
tal angular momentum A^co2 > 0) with experiment!^ and with recent electronic coupled 
cluster /vibrational configuration interaction calculations.'^^'^ Each eigenstate 



Vs,vl,Va) is 

labeled using the quantum numbers of the symmetric stretch Vs, the bend v^ (with I indi- 
cating the vibrational angular momentum, A''co2 > 0; ^^^ the antisymmetric stretch Va- 
The calculated fundamental frequencies of the infrared active bend {uf, = 668.6 cm~^) and 
antisymmetric stretch {ua = 2350.6 cm^^) are accurate to within 1.5 cm~^. The zeroth or- 
der symmetric stretch frequency^ cu^ ~ 1333 cm~^ is about twice as large as the bending 
frequency cob, and the two modes are involved in the accidental Fermi resonance.'^ As a 
result, the vibrational spectrum is organized in polyads with the polyad quantum number 
P = 2fs + Vb] states with P = 2, 3 and 4 are given in Table |TT} In the original version of the 
PES, called 'PESl' in Table|ll| the energies of states (1, 0°, 0) and (0, 2°, 0), belonging to the 
lowest polyad P = 2, are underestimated by 20cm~^ and the difference with the observed 
energies grows rapidly with P. This systematic discrepancy is substantially diminished by 
shghtly rescaling the symmetric stretch, i?+ = (_Ri + i?2)/v^, and the bend aoco via 



R+ -^ V2Rie + {R+ - V2Rie) ■ 1.023 
aoco ^ 180° + (aoco - 180°) ■ 1.0035 . 

The vibrational energies in the scaled 'PES2' agree with their experimental counterparts to 
within 7cm~^ and for the most states below 3000 cm~^ the accuracy is better than 3cm^^. 
The results outperform even the highly accurate calculations of Refs. \TE\ and \W\ shown in the 
second column of Table 11, making 'PES2' one of the best available ab initio potentials of the 
X^S+ state. Since the coordinate dependent dipole moment fix has also been calculated, 
the ab initio intensities of the infrared rovibrational transitions can be directly evaluated.'^Zl 
The other two isomers in Table [T] have never been detected in the gas phase, and the only 
reference data stem from the previous ab initio studies. For the bent OCO, discovered by 
Xantheas and Ruedenberg,'^ the present calculations confirm the C2V symmetric equilibrium 
with the CO bond lengths of 2.51 oq and the OCO bond angle of 73.2°. This minimum is 
located 6.03 eV above the global one, again in good agreement with the previous findings .'^^'^ 
The fundamental excitations in the OCO well, calculated for A^co2 = 0, are Ua = 680 cm~^, 
Ub = 720 cm^^, and Ug = 1550 cm~^. For the linear COO, the calculated CO and 00 
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bond lengths are identical to the ones given in Ref. [2H both are elongated compared to free 
diatoms (2.2 oq vs. 2.14 oq for CO and 2.45 Oq vs 2.28 ao for 00). The calculations place the 
COO minimum at 7.35 eV, about 0.1 eV below the lowest dissociation threshold. 
The ground electronic state correlates adiabatically with two dissociation channels, 

CO2 + /iw(^ph. > 7.41 eV) -^ OCD) + CO{X^i:+) (1) 

CO2 + huj{Eph, > 11.52eV) ^ CC^P) + OalX^S") , (2) 

and the ZPE corrected dissociation energies Dq are shown in Table [ll In channel ([I]), the 
calculated Dq is 0.13 eV less than the experimental value. The deviation might reflect a 
large basis set superposition error introduced by the diffuse functions and as such is the 
downside of the highly accurate vibrational spectrum in Table II The error is independent 
of the arrangement channel, and Dq in channel ^ [closed between 120 nm and 160 nm] is 
equally underestimated. The calculations of Hwang and MebeP^, using a noticeably smaller 
basis set, perfectly agree with the experimental dissociation energy for this channel. 

One-dimensional (ID) cuts through the ground state PES are given for several aoco 
angles in panels (a,c,e) of Figs. [I] and [2] Black solid circles are the raw adiabatic energies. 
The O + CO limit is reached smoothly and no barrier is detected towards the asymptote for 
any orientation of the CO diatom. The same is true for the C + O2 channel, as illustrated 
in Fig. W[e); the potential well in Fig. W[e) is the COO isomer. Angular dependence of the 
X^A' state is shown in panels (a,c,e) of Figs, p^ and k^ for two sets of flxed CO bonds. In 
Fig. |3| Ri is flxed at the FC value; in Fig. |4| it is flxed close to the equilibrium of the bent 
OCO. Consequently, although the carbene-type minimum is perceptible in all panels, it is 
best seen in Fig. kj As CO2 bends, the adiabatic X^A' state (black dots) forms a sharp 
narrowly avoided crossing with the state 2^A' around aoco = 100° (see, for example. Fig. 
[3]). A dynamically meaningful representation for such nearly degenerate pairs is diabatic 
rather than adiabatic. In fact, the black line in Figs, [ij — 14| depicts the X^A' state locally 
diabatized at bent geometries as described in Sect. |IVA[ This is the reason why the black 
dots sometimes switch away from the black line and why lines of different colors cross. In this 
locally diabatic picture, the bent OCO minimum correlates with the state 2M' [purple line]. 
The transition state separating the bent and the linear minima, analyzed by Xantheas and 
Ruedenberg^S and Hwang and MebeP^ is thus the signature of this two-state intersection. 
The two-dimensional (2D) contour plots of the PES of the locally diabatic X^A' state, used 



in the above calculations of vibrational states, are shown in Fig. [5j 

B. Overview of the excited electronic states 

Potentials of the excited electronic states are shown in Figs, [l] and |2] along one CO 
bond for A' and A" symmetries. As with the X state, the solid circles indicate ab initio 
adiabatic energies. In the CO + O arrangement channel, which is the focus of the present 
investigation, all calculated states but one converge to the dissociation threshold ([I]). The 
state 2M", correlating with 1^S~ at the FC point, reaches the higher lying threshold 

CO2 + /iw(Eph. > 11.46eV)^0(^P) + CO(a=^n). (3) 



In the C + O2 arrangement channel [Fig. ^e,f)], three electronic states, the X state and 
the two components of the 11 state, correlate with channel pi). Three other states (E~ and 
A) converge to the electronically excited fragments: 

CO2 + huj{Epi,, > 13.75 eV) -^ CCD) + O2CA) . (4) 

Topographic hallmarks of the excited states can be exemplified using the states 2, 3M'. 
In the FC region near linearity [Fig. U[a)], the states 2M' and 3M' form two sharp avoided 
crossings near Ri = 2.2 oq and -Ri = 2.8 oq. These crossings are in fact two CIs between states 
l^Iig and l^At,.'^'^'^'^ The CIs are not independent: They are connected into a whole line, 
called a 'CI seam' ,1221 unusual properties of whichPare discussed in Sect. IIIC As the molecule 



bends, the gap between the adiabatic states grows. In the lower state 2^A', the intersection 
cone first turns into a broad barrier along the dissociation path [aoco > 170°, Fig. [IJ^c)]. As 
ttoco decreases further, the local minimum near 2.3 ao deepens and the dissociation barrier 
disappears [Fig. [l](e) and Fig. [2]^a)]. In linear COO [aoco = 0°, Fig. |2]^e)], the sharp avoided 
crossing between states 2, 3M' reappears again, this time at Ri = 3.5 oq. Similar to X, the 
state 2M' supports a COO intermediate, although the local minimum lies 0.7eV above the 
C + O2 threshold pi). The evolution of the uppermost state 3^A' with decreasing angle is 
different, because its topography is very much influenced by a pronounced barrier located 
outside the FC zone near Ri = 3.6 qq and separating the fiat inner region from a steep 
decline towards the asymptotic limit. This barrier is distinct over a broad angular range 
and its sharpness suggests a CI with a higher lying state. The nature of this intersection 



becomes apparent in Sect. HIE discussing valence/Rydberg crossings. 



The A" states are similar in many respects. Near linearity, both symmetries mirror the 
topography of the orbitally degenerate states l^Hg and 1^A„, the states 1, 3^A" are involved 
in the same CIs in the FC region [Fig. [lib)], and their behavior in the bent molecule up 
to aoco = 120° closely follows that of 2,3^A'. The l^A" state stabilizes with decreasing 
ttoco, while the 3M" state features a pronounced barrier around -Ri = 3.6 oq [Fig. [l[d,f) 
and Fig. |2](b)]. Similarities persist to aoco = 0° [Fig. [2[f)]: The states 1,3M" form a CI 
at Ri = 3.5 Oo, and the state iM" supports a local COO minimum. Differences with the 
A' symmetry are due to the state 2^A", which in the FC region correlates with l^S" and 
which has no counterpart among the calculated A' states. This state, which at linearity is 
accidentally degenerate with I'^Au, has a single minimum near 2.4 ao and approaches the 
dissociation limit ^ without a barrier. This simple shape is preserved through most cuts 
in Figs, ^andp] The state 2^A" also supports a COO minimum [Fig. ^f)]. 

Potential cuts along bending angle are shown in Figs. |3] and |4} The doubly degenerate 
l^IIg and l^Au states split into A' and A" components for aoco < 180°. Evolution of the 
electronic energies with decreasing aoco was analyzed by Spielfiedel et al. using Walsh 
rules.'^ Based on this analysis, the adiabatic excited states are commonly classified as bent 
or linear. The energy of states 2^ A' and l^A" lowers as the molecule bends, while the energy 
of states 3^ A' and 3^ A" grows [see Fig. ^a,b)]. The global minima of the 'bent' states lie 
near 120° (2^ A') and 130° (iM") and are located below the dissociation limit Q [Fig.|3|a,b) 
and Fig. |4[a,b)]. Bent equilibrium of the state 2^ A', which at C2v geometries becomes ^-82, 
was predicted by Dixon in his analysis of the CO flame emission bands.l^ Finally, the state 
2^ A", which together with 3^ A' and 3^ A" is 'linear', is the least anisotropic of all states in 
a broad vicinity of aoco ~ 180°. 

Properties of various stationary points in the PESs of excited electronic states are given 



in Table III The C^ point group notation is used to label adiabatic states; -Doo/i labels refer 
to the diabatic states at the FC point. Experimental reference data for the excited elec- 
tronic states are scarce and fit into the footnote a. For the bent states, the experimentaP^'^ 
equilibrium CO distances and the bending angle are reproduced within 0.1 qq and 3°, respec- 
tively; the bending frequencies are accurate to within 30cm~^ or better, and the calculated 
band origins lie within the experimental uncertainties. Results of the previous ab initio 
studies^l^ are also shown in Table III Agreement in the equilibrium geometries is excellent, 
with the exception of the intersection-ridden uppermost states 3^ A' and 3^A", for which 
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C2i,-restricted calculations of Refs. [9] and [10] miss the local Cg minima in the upper CI 
cones. Vertical excitation energies T„ are close to those of Ref. [TUl Slight differences are 
not surprising for near degenerate states whose ordering is sensitive to the details of the 
ab initio set up. Peculiarity of the spectral region 120 nm — 160 nm is that T^ values are 
poor approximations to the positions of the absorption maxima because all transitions are 
electronically forbidden and the TDMs at the FC point are strictly zero. Finally, the ab 
initio dissociation energies Dq in channels ([I|, ([2|, (|3|, and Q, also given in Table III, agree 



with the experimental values^'^ within ~ 0.15 eV irrespective of the particular arrangement 
or electronic channel. 

C. Topography of state intersections in the FC zone 

Intersections of electronic states in linear CO2 follow simple 'symmetry rules ^^^'^ which 
severely restrict the intersection topography. Two types of intersections with regard to their 
symmetry properties are prominent in Figs. flta,b) andlslajb): 

(1) Renner- Teller (RT) glancing intersections^ED involve the A' and A" components of the 
l^IIg or the l^Au state. In the rotating molecule, the A /A' interaction ~ Afi/sin^aoco 
is proportional to the projections A and Vt of the electronic (L) and the total angular 
momentum (J = iVco2 + L) on the molecular axis and diverges for aoco ~^ 180° .'^'^^^'^ The 
l^IIg ^ -^^'^t transition is best classified as linear-linear, and the 1^A„ ^ -^'^'^t transition 
is linear-bent. The 'degeneracy manifold' for the intersecting states is the whole (i?i,i?2) 
plane defined by the condition aoco = 180°. 

(2) Two nested CIs involve the 2, 3M' and 1, 2M" or 2, 3M" states and stem from the acci- 
dental ^Ilg/^ A^ crossing. Both degeneracies are lifted linearly along the tuning and coupling 
modes spanning the common branching space.— According to Fig. Q[a,b), the tuning mode 
is -Rco, i-e. a combination of the symmetric and antisymmetric stretch (their irreps are crj"„ 
in the Dooh group); aoco breaking the linear symmetry (irrep tTu) is the coupling mode. As 
a consequence of the 'symmetry rules', the CIs occur along a line Fci{R\,R2) = in the 
{Ri, R2) plane at aoco = 180°: The 'degeneracy manifold' is a ID seam.HESEl 

The CI seam, constructed separately for the 2,3M' and 1,3M" states on a fine {Ri,R2) 
grid,'^ is depicted in Fig. [6J It has two remarkable properties. First, the intersection along 
the seam is fivefold. Two CIs and two RT intersections imply four degenerate states. The 
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hitherto ignored state 1^S~ closely follows l^A^: The A/E energy gap falls consistently 
below 300 cm"^ which is at the limit of the ab initio accuracy. Thus, the l^IIg/l^ A^ and the 
l^Au/l^S^ pairs cross at the same CO bond distances and the total degeneracy is five. Sec- 
ond, the calculated seam traces out a closed loop. Closed CI seams are rarely encountered,'^'^ 
although arguments have been devised to prove their ubiquity.'^^'^ A technical implication 
of closed or strongly curved seams is that local diabatization schemes faiP and global or 
semiglobal diabatizatiorP^ becomes necessary. The impact of closed seams on photodisso- 
ciation dynamics has never been systematically investigated.'^ In CO2, only a small portion 
of the loop is directly accessible to UV light. Nevertheless, paper II demonstrates that the 
seam topology deeply affects the observed absorption spectrum. 

The curved intersection seam is responsible for a peculiar shape of the potentials in the 
FC region. The adiabatic PESs of the states 2, 3M' in linear CO2 are shown in Fig. |6| 
The lower adiabatic state 2^A' has a C2V minimum at Ri = R2 = 2.41 oq, and two Cg 
symmetric saddles near Ri ^ 2.8 ao or R2 ~ 2.8 ao, separating the C2V minimum from the 
dissociation asymptotics. Outside the area enclosed by the seam line, the 2^A' state has 
the n, inside the A character. The saddles hide portions of cusp lines, along which the 
states intersect and which are washed out by the spline interpolation. The topography of 
the upper adiabatic 3^A' state is a literal mirror image of the lower state; one finds a C2V 
saddle at Ri = R2 = 2.41 a^ and two Cg symmetric minima near [Ri, R2) = (2.84 oq, 2.27 Oq) 
and (2.27 ao, 2.84 ao). The state character changes from A outside the seam line to 11 inside. 
The two minima are cusp-like and correspond to the upper cones of the CIs. The two other 
saddle points in the 3^A' state, lying outside the FC zone near (i?i, i?2) = (3.6 ao, 2.2 ao) and 
(2.2ao, 3.6ao), result from an avoided crossing with a higher lying state. 

Changes of the electronic character of A" states across the intersection can be moni- 
tored using matrix elements of the electronic angular momentum L^r^The matrix elements 
{iA"\L;^\jA') along the line passing twice through the closed seam are depicted in Fig. 7 
for z = 1,2,3 and j = 2,3. The states |2M') and \3^A') act as 'probes' whose assignment in 
terms of 11 or A is known. Outside the seam area, most matrix elements are close to integer 



values of 0, 1, and 2. For Ri < 2.2 ao for example. 



{1A"\L,\2A') =land {2A"\L,\2A') 







so that \1A") is a n state, while \2A") is a S state. {3A"\L^\2A') 



vanishes too, but the 



state \3A") is a A state, as confirmed by the matrix element {3A"\Lz\3A') 



2. As the 



CI seam is crossed at Ri = 2.25 ao into the area interior to the seam loop, the three-state 
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intersection induces violent changes in the electronic labels of the adiabatic states. The 
state \3A") becomes a 11 state ( (3y4"|I/2|3A') = 1); the state \2A") acquires A character 
( {2A"\L,\2A') ^ 2), while the state \1A") acquires S character ( {1A"\L^\2A') ^ 0.5), and 
the non-integer values of the latter two reflect strong mixing of the near degenerate A^/E" 
pair. The next reshuffling occurs as the seam loop is crossed outwards at Ri ^ 2.8 Oq. The 
state \1A") becomes 11 again, while the matrix elements involving \2A") and \3A") vary until 
Ri ~ 3.6 Oo is reached, where {3A"\Lz\3A') finally vanishes indicating that \3A") emerges 
from the crossing region as a S state converging towards the upper threshold (|3]), while 
\2A") becomes A state correlating with threshold ([I|. 

The positions of cusps, minima, and saddles in the PESs of the A" states at linearity are 
identical to those in A' states. Exception is the PES of the state 2^A", correlating in the 
PC region with 1^S~. Its appearance, illustrated in Pig. ^ is very much simplified by it 
being tightly linked with the diabatic l^A^ state. The PES has a single C2v minimum at 
Ri = R2 = 2.41 ao, is free from cusps seen in other adiabatic PESs, and provides a nature's 
illustration of the shape of the diabatic 1^A„ PES. 

CIs at linearity can also be recognized in the potentials plotted in the (i?i, aoco) plane. 
An example including four electronic states at R2 = 2.2 oq is shown in Fig. |9j Although 
the characteristic features are no longer C2V symmetric, a maximum in the 2^A' state at 
Q^oco = 180° and -Ri = 2.8 qq and two minima in the 3^A' state at aoco = 180° and 
Ri = 2.3 ao and Ri = 2.8 oq are recognizable in panels (a) and (b). Close to linearity, the 
A" states [panels (c) and (d)] maintain the same contour maps as A' states. 

D. Intersections and avoided crossings in bent CO2 

As CO2 bends and all above degeneracies are removed, the primary topographic features 
(cusps, minima, and saddles) remain clearly visible up to aoco ~ 160°. This is demonstrated 



in the (i?i,aoco) plane in Fig. [9] and in the {Ri, R2) plane in Fig. 10 Further decrease 
of aoco spawns new avoided crossings, a detailed map of which is given in the angular 
cuts in Figs. |3] and |4| All pairs of calculated states become near degenerate at various 
geometries. For example, the X state forms avoided crossings successively with states 2M' 
at aoco = 100° and 3^A' at aoco = 80° (Fig. III). The state 3M" approaches closely 2^A" at 



aoco ~ 120° - 130° and VA" at aoco = 90°, and the states VA" and 2M" cross between 70° 
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and 80°. Solid lines in Figs. [3] and |4] cross because they refer to the states locally diabatized 



through the 'bent' avoided crossings (see Sect. IV A). Potential curves of the states 3^A' and 
3^A" bear evidence of strong interactions with the next higher states. Especially pronounced 
in Fig. p[d,f) and Fig. p[b,d) are the sharp near-intersections in 3^A" around 155°. 

As a result of multiple avoided crossings, local bent equilibria are found in all calculated 
states. The diabatic origin of a particular local minimum is invariably different from the 
adiabatic one. Consider the carbene-like bent OCO with aoco ~ 70° in Fig. |4[a). Purely 
adiabatically (solid dots), this minimum belongs to the ground electronic state. In the 
locally diabatic picture, bent OCO belongs the state 2M' (purple line). An avoided crossing 
between states 2^A' and 3^A', recognizable in panel (a) around aoco ~ 95°, implies that 
another diabatization might re-assign the bent OCO minimum to the 3^A' state (brown 
line). Furthermore, the broad barrier in 3M' near aoco ~ 140° indicates another avoided 
crossing with the next higher A' state which thus would be the true owner of the bent OCO 
minimum. Similar analysis applies to the bent OCO in the A" states in Fig.llFb). Depending 
on the chosen representation, it can be ascribed to the fully adiabatic state iM" (dots), to 
the locally diabatic state 3^A' (brown line), or — via the sharp near- intersection around 155° 
— to the next higher A" state. Other bent conformations in the excited states result from 
state interactions, too. One such isomer with the valence angle of 100° is created via an 
avoided crossing in 3^A' [Fig. p[c,e)]. Crossing of the same 3^A' state with the next higher 



state between 60° and 70° leads to another high-energy bent minimum [Fig. |4Fa,e)]. 

E. Beyond the first six states: Valence/Rydberg crossings 

Many local barriers and minima in the states 3^A' and 3^A" are due to intersections with 
'invisible' higher lying electronic states. Previous detailed ab initio studies exposed these 
'invisible' intersection partners as mainly Rydberg states.'^'^^'^^ The aim of this section is 
to illustrate how the Rydberg/valence interaction affects the potential profiles along the 
dissociation path and to outline the geometries at which valence states can effectively drain 
population from the Rydberg manifold. To this end, CASSCF calculations of the first 10 
states, 1 — 5^A' and 1 — 5M", have been performed with the d-aug-cc-pVQZ basis set for the 
CO bond distances R2 G [1.7ao, 5.0ao] and the OCO bond angles aoco ^ [100°, 179°]. One 
CO bond is kept fixed at Ri = 2.4 ao, so that only CO + O arrangement channel is covered. 
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The CASSCF potentials along i?2 are shown in Fig. 11 All five A" states and four A' 



states converge to thresholds ([I| or (Is]). The states correlating with the lowest threshold ([I]) 
are of S"*", A, and 11 symmetry. The states correlating with the highest threshold pi) are 
S^, A, and 11. One A' state converges towards the asymptote 

CO2 + huj{Epi,, > 9.64eV) -^ 0(^5) + CO(X^S+) . (5) 

At linearity, this state is the Rydberg state l^S+(37ru).'^ Compared to Fig. [l| many new 
avoided crossings are found in Fig. [ll|(a,b). The Dooh notation of a state \i) in these panels 
is related to the expectation value {i\Ll\i) which a diabatic state preserves across the inter- 
sections. Potential curves in Fig. [ll[a,b) are color coded according to the (L^) values: S 
states {{Li) = 0) are shown in green; 11 states {{L^) = 1) are blue, and A states {{L^) = 4) 
are red/orange. Clearly, most adiabatic curves change color more than once as the CO bond 
stretches: They are tailored out of several diabatic states. 



Familiar from the discussion in Sect. IIIC are the fivefold crossings involving 1^11^(50" 



gj 



and the accidentally degenerate pair l^A„(27r„)/l^S^(27ru). These crossings are marked with 
arrows in Fig. [ll[b). As in the MRD-CI calculations, the A„/E~ pair is easily recognizable 
using matrix elements of L^. At linearity, (L^) for these two states is either or 4, but 
already a tiny deviation of 2° causes {L^) to collapse towards an average value of ~ 2. 
This state mixing is stressed in Fig. [ll[a,b) with the A^j/S" label and with the red/green 
color. The second edition of the mixed A^/S^ pair is found around 11.5 eV. The Rydberg 
state 2^Au(37r„), shown with an orange line in Fig. [ll[a,b), carries as a satellite the state 
2^S^(37ru) [green line in Fig.[ll|(b)]. Again, (L^) for these states assumes a non-integer value 
between 1.5 and 2.5 in even slightly bent molecule. 

The valence pair lA^/lS" converges towards the uppermost threshold ^ and succes- 
sively traverses the higher lying states. The prominent dissociation barrier in the 3^A' and 
3^A" states near R2 = 3.6 oq is a clear signature of the sixfold valence/Rydberg crossing 



involving 1A„/1S~, the l^S^ state [green line in Fig. 11 ^a)] and a repulsive -R^A„ state 
[orange line in Fig. [ll[a,b)]. l^S^ is the optically bright Rydberg state responsible for the 
strong absorption band around 11.1 eV (lll.Tnm).'^ The state R^A^ corresponds to a pair 
of strongly repulsive M'/M" states descending towards threshold ([T]) from very high ener- 
gies. Two sections of this repulsive potential curve are seen in Fig. [lT|^a,b) in the intervals 



[3.2ao, 3.7ao] and [3.7ao, 5.0ao]. Thus, the asymptotic repulsive portions of the 3 A' and 
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3M" states diabatically belong to -R^A„. The sixfold crossing lS+/_R^A,u/lA„/lS~ near 
i?2 = 3.6 ao is not the only one involving -R^A„ state. At higher energies and shorter CO 
bonds, R^Au intersects the Rydberg pair 2A„/2S^ (cf. Fig. [ll[a,b) near R2 ~ 3.0 Oq). 
The diabatic -R^A^ state distinctly stands out because it strictly preserves the projection 
{LD ^ 4 even for aoco = 165°, while (L^) values for the other states become non- integer. 
The above discussion is valid for both A' and A" states. The A" symmetry block in Fig. 



11 b) contains one more state, namely the Rydberg state 2^Ilg{6ag) missing among the A' 
states where it would have been 6A'. This state, materializing out of nowhere at R2 = 3.0 og, 
is involved in the i?^Au/2Au/2S^ crossing — making it a sevenfold intersection. 

While the adiabatic gaps in the valence/valence intersection in^/lAu/lS^ grow as aoco 
deviates from 180°, all Rydberg/Rydberg and Rydberg/valence intersections not only remain 
recognizable in bent CO2, but clearly sharpen up. This is illustrated in panels (c) — (f) 



of Fig. 11 drawn for aoco = 175° and 160°. As a side result, the repulsive R^A' and R^A" 
states, deriving from R^A^, remain distinct at all angles, despite new 'bent' intersections in 
the 3,4,5^^4' and 3,4, 5^y4" states [see, for example. Fig. [ll|(e,f)].^ Via these intersections, 
carbon dioxide excited with vacuum UV light can — at bent geometries — reach any of 
the shown dissociation channels. Although a detailed analysis of these multiple pathways 
is beyond the scope of the present work, a brief description of intersections involving the 
optically bright ^S+ state is added at the end of this section. 

Potential curves along aoco are shown in Fig. [12} The crossing patterns, expecially 
numerous among A" states, explain intricate topography of states 3^A' and 3^A" in Fig.|3} For 
example, the sharp barrier in the state 3^A" at angles between 150° and 170° [Fig. p[d,f ) and 
Fig. [3tb,d)] originates from a complicated three state crossing around 165° involving states 
3,4, 5M" [Fig. [l2|(b)]. Diabatically, the decreasing branch of the 3^A" state at aoco < 160° 
belongs to 4M" — the state which at linearity merges into the Rydberg pair 2Au/2S~. 
This diabatic state can be followed to even smaller angles through another intersection, this 
time with 2^A" near 130°. Ultimately, as indicated with an arrow in Fig. |l2|(b), it is this 
diabatic state which the carbene-like OCO minimum belongs to. The avoided crossings in 
the A' symmetry block are similar and simpler. The three state crossing of 3,4, 5M" is a 
mere ghost because the adiabatic gap is almost 2 eV wide and the barrier in the 3^A' state is 
broad and low. The state 4M', originating from 2^Au, stabilizes upon bending and is easily 
traced to smaller angles through the intersection with 3M' near aoco = 150°. Again, the 
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carbene-like OCO correlates diabatically with the state originating from 2^Au at hnearity. 

The analysis of broadening and splittings in the strong absorption band of the Rydberg 
state 1^S+ focused on the interactions of 1^S+ with valence states.'^'^^J^ Present calculations 
reveal numerous avoided crossings with both valence and Rydberg states inside and outside 
the FC zone. At linearity, the 1^S+ state cuts twice through the Rydberg state 2'^Au [Fig. 



11 'a)], and both crossings persist to smaller angles [Fig. ll|^c,e)]. Another crossing occurs 
near R2 = 3.7 Oq (for all calculated angles) and mixes 1^S+ with -R^A„. This interaction 
strongly perturbs the ab initio (L^) values but vanishes beyond Ri = 4.0 oq. The potential 
cut along aoco in Fig- [l2|a) exposes another avoided crossing at ~ 120° which involves the 
states 4, SM' and leads to a local minimum in SM' (the state, correlating with l^S^ at these 
Rco)- Finally, the barrier in 5M' near 150° implies interaction with the next higher state, 
which according to the analysis of Ref. H2] correlates with 2^S+ at linearity. While all these 
crossings can redirect population from the optically bright l^S^ state along various linear 



and bent routes. Fig. 11 clearly demonstrates that the diabatic l^S^ state, calculated at 



the CASSCF level, is repulsive at all geometries and thus can dissociate directly. 

IV. QUASI-DIABATIZATION OF THE VALENCE STATES 

As explained in the Introduction, the diabatic representation,-^- although not generally 
indispensable, is best suited for the particular implementation of nuclear quantum dynamics 
used in paper II. Rigorously speaking, diabatization has to be performed simultaneously on 
all six calculated valence states, because each pair of states intersects either in the FC zone 
or at bent geometries. Simplifications to this 'Herculean taskSS stem from the expectation 
that the two groups of intersections influence photodissociation in different ways: While 
the 'bent' CIs affect later stages of the product formation, the electronic branching ratios, 
and/or the rovibrational photofragment distributions,'^ the 'FC CIs are directly responsi- 
ble for the shape of the observed absorption spectrum. This distinction guides the practical 
construction of the diabatic states: The 'bent' CIs and the 'FC CIs are analyzed at two 
different levels of detail, commensurate with their expected impact on the absorption spec- 
trum. As a result, the complete multistate problem splits into several steps and the need 
for a global diabatization of six multiply intersecting states is obviated. 
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A. Intersections at bent geometries: Local diabatization 

The vicinity of 'bent' CIs is diabatized locally, using the energy-based scheme as described, 
for example, by Koppel.!^ All considered CIs are shown in Figs. ^ and E^ and include 



(a) the XM72M' pair at aoco = 90° - 110°; 

(b) the 2M'73M" pair at aoco = 120° - 150°; 

(c) the IM'VSM" pair, with 3M" state diabatized in step (b), at aoco = 70° - 100°. 

The XM'/2M' pair is considered as an example. At C2V geometries, these states belong 
to Ai and B2 irreps, so that their CI is 'symmetry allowed'. 1^ Its branching space includes 
Q^oco as the tuning mode (irrep Ai) and the antisymmetric stretch R_ as the symmetry 
breaking coupling mode (irrep i?2)- The 2x2 diabatic potential matrix V^ is constructed 
on the ab initio grid of bond distances (-Ri, R2) from the diagonal adiabatic potential matrix 
V" using the orthogonal adiabatic-to-diabatic transformation (ADT/^ 

yd _ s^V^S , (6) 

with the transformation matrix S parameterized by the mixing angle 0, 

r>i^D D ^ -*- ^- /ttoco -"oco(-Ri'-^2)A /-x 

e(i?„i?„aoco) = -2 arctan ^ WoiR^.R^) ) ' ^^^ 

which varies between and 7r/2; aQ^^Q^Ri, R2) denotes the crossing point of the diabatic 
potentials (9 = 7r/4); the width Wq of the function 0(aoco) is evaluated from the adiabatic 
energy gaps AV^/2 and the average slope of the diabatic potentials F, 



Wo 



2F 



(8) 



at each grid point {Ri, R2) for agco — 100°. With the off-diagonal coupling in Eq. ^ 
neglected, the ground state X'^A' becomes completely decoupled from other A' states. The 
remaining diagonal matrix element in V^ is the potential waiting to be further diabatized 
at linearity. The diabatization is kept local by restricting the interval, over which 0(aoco) 
varies from to tt/2, to the vicinity of aQco- 

e(i?i,i?2,aoco) =5-1 [5'20(^i,^2,ttoco) + 2(1 ~5'2)J , (9) 
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The two switching functions 5'i,2(«oco), 

fl'i,2(aoco) = [1 + exp(-(aoco - tti,2)/Ai,2)]"^ 



(10) 



with «! = 95°, «2 = 105°, and Ai = — A2 = 5°, are adjusted to restrict the diabatization to 
angles 90° — 110°. If (ai — 0:2) — )■ in (71 2, the so-called 'diabatization by eye' is recovered, 
which corresponds to a relabeling of ab initio energies at the crossing angle Ooco-^^^ 

The main advantage of local diabatization, the results of which are shown in ID in Figs. 
|3] and |4| is that smooth potentials are created for strongly bent geometries at relatively low 
cost. The quality of the resulting PESs in 2D, either in the {Ri,R2) or in the (i?i,aoco) 
plane is illustrated for all states in Figs. [5} [8} |9| [10} and 13 For the lowest excited states in 
each symmetry block, 2^A' and iM", the global bent equilibrium at 120° or 130° [Fig. ^a,c)] 
is not affected by the local diabatization. In contrast, the carbene-like bent OCO minima 
change the owner: They appear in the state 2M' [local minimum near 70° barely visible in 



Fig. |9[a) but pronounced in Fig. 13;^a)] and in the state 3^A" [Fig. [l3|(d)]. Barriers and local 
minima arising in the avoided crossings with Rydberg states are not diabatized within this 
scheme. Examples are the sharp barrier in 3M" near 150° [Fig. ^d) and Fig. [l3[d)] and the 
deep local minimum in 3^A' near 100° [Fig. ^b) and Fig. Q3[b)]. 



B. Intersections at linear geometries: Regularized diabatic states 

In the second step, the CIs in the FC region are diabatized; kinematic singularities due to 
glancing intersections between A' and A" states are not considered. Accurate diabatization 
schemes rely on direct numerical differentiation of the electronic wave functions with respect 
to nuclear coordinates or on the orbital rotation methocP^'^as implemented in MOLPRO.'^ 
Both types of calculations become prohibitively expensive with the d-aug-cc-pVQZ basis 
set, and further approximations are needed in order to find the diabatic potential matrix, 
consisting of a 2 x 2 block of A' and a 3 x 3 block of A" states: 



V° 



( Va' 


^n'A' 






\ 


^n'A' 


Vw 












^S" 


Vs"n" 


Vs"A" 






Vj]"u" 


Vu" 


Vn"A" 


\ 




Vs"A" 


Vn"A" 


Va" ) 



(11) 
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The five diagonal matrix elements are the five diabatic states S", IT', 11", A', and A", 
with the smooth intersecting potentials in the {Ri, R2\aoco = 180°) plane constructed to 
coincide at aoco = 180° with the ab initio PESs for l^S^, l^Hg, and l^A^,, respectively. 
The semi-global diabatization scheme, akin to the vibronic coupling modeP^'^and adjusted 
to the topography of the closed fivefold CI seam, has already been introduced in Ref. |3l 



Diabatization proceeds in two steps. First, a model diabatic matrix of the form Eq. (11) 
with elements V^j is constructed. Due to orbital degeneracy, Vy^, {Ri, R2) = V^„ (-Ri,i?2) 
and V^, {Ri, R2) = V^J{Ri, R2). Moreover, the accidental near-degeneracy of states 1^S~ 
and 1^A„ implies Vj^„ {Ri, R2) ~ V^„ (-Ri, -R2) in a broad vicinity of the closed CI seam — the 
property which substantially simplifies diabatization of A" states. Deviations from linearity, 
measured by the coordinate Qu ~ sin aoco, are included in the model via off-diagonal matrix 
elements represented as symmetry adapted expansions in Qu- 

N' 

vSl = Y.^k{Ri,R2)Ql'^' (12) 

fc=0 

N" 

^E"n" = ^n"A" = ^/^fc(-Ri,-R2)Q„ "^ • (13) 

fc=0 

Couplings of the accidentally degenerate S" and A" states to 11" are set equal, while the 
matrix element V^//^// for the RT-like S"/A" interaction is neglected. The model is complete 
after the expansion coefficients in V^ ■ are calculated from a non-linear least-squares fit to ab 
initio energies. In the second step, the regularized diabatic states approach is invoked, ^^ and 
the matrix elements V^, are used to define the orthogonal ADT matrix, which is applied to 
the adiabatic matrix V" via Eq. ([6]) giving the desired diabatic matrix elements of Eq. (11) 
on the full ab initio grid. Final interpolation is performed using 3D splines. 

The diabatization is localized to the vicinity of the CI seam by modifying matrix elements 



V^A°)(*^j)mEq. (|12| 



y-.(°) = 'A (14) 

'^ l+exp[(p-po)/Ao]' ^ ^ 

with p = [{Ri — 2.605 flo)^ + {R2 ~ 2.605 cto)^] ^ being a radial distance from the center of a 
circle enclosing the seam; po = 1-0 oq and Aq = 0.1. Adiabatic and diabatic states are forced 
to coincide if either bond becomes longer than ~ 3.8 ao, so that diabatization regions for 
linear and bent CIs are cleanly kept apart. 
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The constructed representation describes best the vicinity of hnearity in which non- 
adiabatic transitions occur. The exact range of vahdity is determined by the length of the 



expansion in Eq. (12). The choice A^' = A^" = 1 gives a model which fits ab initio data with 
a root mean square error of ~ 180 cm~^ for angles aoco > 150°. The ADT, constructed 
using this model, is guaranteed to remove kinematic singularities at the CIs, but leaves the 
strength of residual non-adiabatic couplings unspecified.'^ The ultimate test of the scheme 
is the quantum mechanical absorption spectrum described in paper II. In order to assess 
the accuracy of the truncated Qu expansion, three diabatic representations are constructed, 
based on the expansion coefficients obtained from fitting in three different angular ranges 
180° — 170°, 180° — 160°, and 180° — 150°. The corresponding absorption spectra are virtually 



identical. The spectra are also insensitive to small variations in po and Aq in Eq. (14) — the 



modifications take place too far away from the crossing seam to affect nuclear dynamics. 



Another test of the constructed ADT is given in Fig. 14 , in which the model mixing angle 
G for states 2, 3M' is compared with the ab initio one calculated using MOLPRO with a 
smaller cc-pVQZ atomic basis set. The dependence O(-Ri) on the CO bond length has a 
characteristic bell shape: The closed CI seam is intersected twice giving rise to the ascending 
and the descending branch. The curve 0(-Ri) fiattens out as aoco decreases and CO2 leaves 
the degeneracy plane at 180°. Agreement between the model and the ab initio results is 
satisfying for all angles. A constant shift of 25° applied to the ab initio mixing angle has no 
effect on the strength of non-adiabatic coupling proportional to dQ/dRi. 

One-dimensional cuts through the diabatic PESs (diagonal matrix elements corresponding 



to states l^IIg and l^A^) are shown in Fig. 15 They cross at all geometries and can be 



directly compared with adiabatic states in Fig. [Tj The off-diagonal coupling matrix elements 
are large in the intersection region vanishing off towards the asymptotic channels. Diabatic 



matrix elements are further illustrated in Fig. 16 in the (i?i, R2) plane and in Fig. 17 in the 



(i?i,aoco) plane. In all cuts, the diabatic PESs smoothly depend on internal coordinates. 
In the {Ri,R2) plane, the off-diagonal diabatic coupling stays localized in the inner region. 
In contrast, the coupling along bending coordinate is delocalized across a substantial aoco 
range in the (-Ri,aoco) plane. As a result, the angular shape of the diabatic potentials 
is distorted compared to the adiabatic case: Diabatic potentials along the coupling mode 
Q^oco are close to the average adiabatic potential ^iV^ + VJ"). 
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V. CONCLUSIONS 

This paper describes properties of global PESs of the first six singlet electronic states 
of CO2 constructed from about 5000 symmetry unique ab initio points calculated with the 
d-aug-cc-pVQZ basis set using the MRD-CI method. The main results can be summarized 
as follows: 

1. Calculations accurately reproduce the known benchmarks for all states and estab- 
lish missing benchmarks for the future calculations: Bond distances and bond angles 
are accurate to within 0.1%, known fundamental frequencies (mainly ground state 
X^A') are accurate to within 1.5 cm~^, the accuracy of the vertical excitation energies 
is expected to be better than 0.05 eV; dissociation energies agree with experimental 
thresholds within 0.15 eV for four covered arrangement channels. 

2. Local equilibria are abundant in the calculated states. Bent OCO isomer is found in 
the adiabatic states 1,3M' and 1,3^A". Linear COO is found in 1,2^A' and 1,2^A". 
Their diabatic electronic origin is clarified, and the properties, including equilibrium 
geometries, excitation energies, and vibrational frequencies, are established. 

3. Near degeneracies can be found for each pair of six valence states, at linear or bent ge- 
ometries, or at both. Avoided crossings and conical and glancing intersections literally 
shape the observed topography of the excited electronic states. Detected intersections 
are not limited to the valence manifold and the search for electronic origins of local 
minima and barriers involves valence/Rydberg and Rydberg/Rydberg intersections. 

4. Characteristic for state intersections in CO2, both conical and glancing, is that they 
include several states. In the FC region, a fivefold intersection between l^S^, l^LI^, and 
1^A„ states is found. The seam of this intersection forms a closed loop, spectroscopic 
manifestations of which are discussed in paper II. Outside the FC region at linearity, 
six- and sevenfold intersections are predicted, some of which persist over extended 
angular range in the bent molecule. 

5. Diabatic 6x6 potential matrix, with all elements smoothly depending on internal co- 
ordinates, is constructed using two-step local diabatizations of linear and bent conical 
intersections. 
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It is tempting to try to infer the course of photodissociation and the principal features of 
the absorption spectrum — the outcome of a comphcated quantum mechanical calculation 
— from the constructed PESs alone. Two issues have to be resolved if one deals with five 
interacting states. The first is the strength of diabatic (intra-symmetry) and RT (inter- 
symmetry) coupling. If the off-diagonal vibronic coupling is weak, the diabatic potentials 
should be chosen as 'zeroth order guides'. If the vibronic coupling is strong, it is the adiabatic 
description which becomes relevant — and the difference between the two pictures is striking, 



especially along the bending angle as Figs. [9] and [17] demonstrate. As discussed in paper 
II, the vibronic coupling is strong, the RT interaction between A' and A" states is to a 
large extent quenched, and the adiabatic potentials can be used for qualitative analysis. 
The global minima of states 2^A' or iM" are bent and lie ~ 2.0 eV below the dissociation 
threshold ([I| or ~ 4.5 eV below the FC point. The route connecting the FC region with these 
bent equilibria is barrierless. In contrast, there is a barrier to dissociation near linearity — 
the leftover of the lower cone of the II/A CI. Thus, one expects the low energy bands in the 
absorption spectrum, associated with 2^A' and iM" states, to reflect highly excited bending 
motion. This interpretation is commonly given in the literatureP Instead of dissociating 
directly, the molecule bends first. With growing photon energy, the contribution of direct 
dissociation through linear geometries will certainly grow, because the barrier is only about 
0.2 eV high and is located ~ 0.4 eV below the FC point. Above ~ 9 eV, the valence states 3M' 
and 3^A" will contribute to the observed spectrum. These 'linear' states are separated from 
the dissociation asymptote by a ~ 1 eV high and broad barrier, with the implication that 
CO2 in these states can decay only through non-adiabatic interactions with the dissociation 
continuum of the lower states. In other words, one expects to see a resonance-dominated 
absorption spectrum. The next qualitative change in the absorption spectrum within the 
valence manifold can be expected after the photon energy reaches the top of the dissociation 
barrier in the upper valence states and allows direct dissociation from linearity. 



The above discussion is based on ID and 2D potential cuts and the data in Table [ITI — 
given the adiabatic representation is the adequate one. However, there is another important 
piece of information still missing, namely the TDMs with the ground state. As has already 



been mentioned in Sect. IIIB the electronic transitions in the wavelength range of 160 nm 
— 120 nm are forbidden. The bands are observed only because the TDMs are not constant 
but strongly change with molecular coordinates as one moves away from the high-symmetry 
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FC point. This dependence is a manifestation of the Herzberg-Teller effectP which plays the 
leading role in shaping the absorption bands, is at least as important as the potential profiles 
discussed above, and has to be considered on equal footing with the state intersections. The 
discussion of the coordinate dependence of the TDMs is deferred to paper II. 
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TABLE I: Properties of the adiabatic ground electronic state X^T,'^: The OCO bond angle (aocOj 
in °), CO and 00 bond distances (in gq), and the energy AE relative to the global minimum (in 
eV, including the vibrational ZPE corrections), calculated for three equilibrium geometries and the 
two arrangement channels. 



Geometry 


"OCO 


-Rco,i 


-Rco,2 


^oo 


AE^ 


Reference 


OCO linear 


180.0 


2.1991 


2.1991 


4.3982 


0.0 


this work 




180.0 


2.2050 


2.2050 


4.4100 


0.0 


m 




180.0 


2.1924 


2.1924 


2.3848 


0.0 


I9I10I 




180.0 


2.1960 


2.1960 


4.3920 


0.0 


m 


OCO bent 


73.2 


2.51 


2.51 


2.99 


5.90 


this work 




72.9 


2.52 


2.52 


2.99 


6.04 


m 




73.1 


2.53 


2.53 


3.01 


5.97^ 


[231 


COO linear 


0.0 


2.20 


4.65 


2.45 


7.35^= 


this work 




0.0 


2.20 


4.65 


2.45 


7.14 


m 




0.0 


2.17 


5.37 


3.20 


6.86 


m\ 


Oi^D)+CO{X^^+) 


— 


2.14 


oo 


oo 


7.28 


this work 






2.175 


oo 


oo 


7.64 


m 






2.132 


oo 


oo 


7.41 


|4|28p 


C(3P)+02(X3S-) 


— 


oo 


oo 


2.30 


11.35 


this work 






oo 


oo 


2.356 


11.49 


m 






CXD 


oo 


2.283 


11.52 


1^ 



"The ab initio ZPEs are: ZPE(C02) =0.314 eV, ZPE(CO) =0.130 eV, and ZPE(02) =0.095 eV. 

''Experimental values. 

"^Energy without ZPE correction. 
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TABLE II: Energies (in cm^^) of the low lying vibrational states in the ground electronic state 
X^S+ measured relative to the ground level (0,0*^,0). The states calculated with the original ab 
initio potential ('PESl') and with the potential rescaled along the symmetric stretch and the bend 
('PES2') are compared with the best theoretical estimates of Refs. [T8)19l and with the experimental 
values of Ref. ITT] (denoted 'Exp'). The energy difference £^observed — -^'calculated (in cm~^) is shown 
in parenthesis. The quantum numbers (vs,vl,Va) are defined in text. 









PESl 


PES2 




State 


Refs.ll8U9| 


Exp 


(0,0°, 0)" 


0.0 


0.0 


0.0 


0.0 


(0,li 


0) 


669.1 (-1.7) 


668.6 (-1.2) 


670.3 (-2.9) 


667.4 


(1,0° 


0) 


1288.9 (-3.5) 


1265.3 (20.1) 


1284.1 (1.3) 


1285.4 


(0,22 


0) 


1339.6 (-4.5) 


1336.6 (-1.5) 


1336.8 (-1.7) 


1335.1 


(0,2° 


0) 


1389.3 (-1.1) 


1373.5 (14.7) 


1389.2 (-1.0) 


1388.2 


(1,1^ 


0) 


1938.0 (-5.5) 


1913.2 (19.3) 


1933.4 (-0.9) 


1932.5 


(0,33 


0) 


2011.4 (-8.2) 


2005.9 (-2.7) 


2006.1 (-2.9) 


2003.2 


(0,31 


0) 


2080.0 (-3.1) 


2061.7 (15.2) 


2079.6 (-2.7) 


2076.9 


(0,0° 


1) 


2349.2 (0.0) 


2350.6 (-1.4) 


2351.4 (-2.2) 


2349.2 


(1,2° 


0) 


2552.0 (-8.6) 


2516.4 (27.0) 


2549.5 (-6.1) 


2543.4 


(2,0° 


0) 


2676.3 (-5.2) 


2626.0 (45.1) 


2668.3 (2.8) 


2671.1 


(0,40 


0) 


2809.1 (-12.0) 


2761.8 (35.3) 


2790.2 (6.9) 


2797.1 


(1,22 


0) 


2589.8 (-4.7) 




2581.6 (3.5) 


2585.1 



"The energy difference between the ground vibrational level and the potential minimum is 2516 cm ^ for 
PESl and 2533 cm^^ for PES2. The ZPE = 2508.5 cm"\ given in the NIST database,^'^ is evaluated from 
the fundamental frequencies via ujb + ^(Ws + i^a)- For PESl and PES2, this value is 2503.6 cm~^ and 
2512.5 cm~^, respectively. 
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TABLE III: Properties of the PESs of the first five excited singlet states of CO2: Vertical excitation 
energy Ty (in eV); band origin Tq (in eV), which includes ZPEs of the ground and the excited 
electronic states; equilibrium CO bond lengths Ri^e and -R2,e (in oo); equilibrium OCO bond angle 
Ue (in °); quantum mechanical vibrational frequencies Ug, uJa, and Uh near equilibrium (in cm^^); 
quantum mechanical dissociation thresholds Do(0/CO) and Z)o(C/02) in the 0+CO and C+O2 
arrangement channels, respectively (in eV). 



Cs 


l^A"" 


2^A"' 


2U" 


31^'^ 


31^// 6 


Dooh 


^n. 


^n. 




iA„ 


^A„ 


0+CO channel 


i^/iS 


1^/1 s 


3p/3n 


1^/1 S 


i^/iS 


C+O2 channeP 


i^/iA 


i^/iA 


iD/iA 


3P/3S 


3P/3S 


To 


5.36 


5.39 


7.95 


8.67 


8.67 


T. 


8.92 


8.92 


8.79 


9.16 


9.17 


Ty, Refs.Iini 


9.00 


9.00 


9.19 


9.28 


9.28 


Rl,e 


2.37 


2.36 


2.41 


2.25 


2.25 


R2,e 


2.37 


2.36 


2.41 


2.80 


2.80 


Re, Rei.mm 


2.38 


2.38 


2.40 


2.29 


2.29 


ae 


127.3 


117.9 


176.0 


180.0 


180.0 


ae, Ref.fTOl 


127.0 


117.8 


180.0 


180.0 


180.0 


UJs 


1283 


1340 


1015 


520 


550 


^a 


905 


865 


1118 


1560 


1550 


idb 


670 


580 


577 


2290 


3100 


^o(0/CO)'^ 


7.27 


7.27 


11.31 


7.27 


7.27 


T'o(C/02)'^ 


11.34 


11.34 


13.59 


13.59 


13.59 



'^Experimental data from Ref. 10 are: For lU" R^ = 2.28 ± 0.02 oq, a^ = 129 ± 1°, Tq < 6.2 eV, v^ ~ 

632cm-i; for 2^' R^ = 2.35 ± 0.015 ao, cxe = 122 ±2°, Tq w 5.7 eV, Vb w 600cm-i. 
''For the i^A' and "i^A" states, the geometries refer to the local minimum in the FC region. The vibrational 

frequencies are strongly perturbed by CI cusps, and the corresponding ZPEs are omitted in Tq- 
"^Dissociation thresholds, labeled with electronic states of atomic/diatomic fragment, are correlated with 

the diabatic states of CO2. 
''Experimental dissociation thresholds are given in Eqs. (nj), ([2|, ([3|, and Q. Ab initio ZPEs are: 

ZPE[C02(XiS+)] = 0.314eV; ZPE[CO(XiE+)] = 0.13GeV; ZPE[C0(a3n)] = 0.099eV; ZPE[02(X3e^)] 

= 0.095 eV; ZPE[02(iAg)] = 0.080 eV. 
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FIG. 1: Cuts through the PESs of the states 1,2,3^' (a,c,e) and 1,2, 3^ A" (b,d,f) along one CO 
bond distance -Rco- The second CO bond is fixed at 2.2 oq; the fixed bond angle is indicated in 
each panel. Dots are the ab initio adiabatic energies. Solid lines are states diabatized at bent CIs 



as described in Sect. IV A For A" states, the ground electronic state X^A' , shown with gray dots, 
sets the vertical energy scale. At linearity, components of the same orbitally degenerate state Hg 
or Au have the same color in panels (a) and (b). 

FIG. 2: Cuts through the PESs of the states 1,2,3^' (a,c,e) and 1,2, 3^ A" (b,d,f) along one CO 
bond distance -Rco- The fixed bond angle is indicated in each panel. The second CO bond is 
fixed at 2.2 ao in (a-d). In (e) and (f), the 00 bond is kept fixed at 2.3 oq. Dots are the ab initio 



adiabatic energies. Solid lines are states diabatized at bent CIs as described in Sect. IV A For A" 
states, the ground electronic state X^A', shown with gray dots, sets the vertical energy scale. 

FIG. 3: Cuts through the PESs of the states 1,2,3^' (a,c,e) and 1,2,3^" (b,d,f) along bond 
angle aoco- The two CO bonds are fixed, one at 2.2 ao and the other as indicated in each panel. 
Dots are the ab initio adiabatic energies. Solid lines are states diabatized at bent CIs as described 
in Sect. HV^ 

FIG. 4: Cuts through the PESs of the states 1,2,3^' (a,c,e) and 1,2,3M" (b,d,f) along bond 
angle aoco- The difference with Fig. [3] is that one CO bond is fixed at 2.5 cq. 

FIG. 5: Contour maps of the ground electronic state X^A': (i?i,i?2) plane, with aoco fixed at 
179° (a) and 70° (b), and (i?i,aoco) plane with Rqo fixed at 2.2 ao (c) and 2.6 oq (d). Energy of 
the dotted contour in (a), (c), and (d) is 2.1 eV, in (b) it is 8.4eV. Contour spacing is 0.35 eV. 

FIG. 6: Contour maps in the (-Ri,i?2) plane of the adiabatic states 2^A' (a) and 3^A' (b) for 
Q^oco = 179°. The states are components of l^H^ and l^A^, respectively. Black solid line in each 
panel indicates the closed CI seam for A' (a) and A" (b) symmetry states. Black dashed line in (a) 



shows the path along which the matrix elements 



{iA"\Lz\jA') in Fig. 7 are calculated. Energy of 



the dotted contour is 7.5 eV, and the contour spacing is 0.25 eV. 
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FIG. 7: MRD-CI matrix elements of L^, {iA"\Lz\jA') , for electronic state pairs (i = 1, j = 2) 
(black), (i = 2,j = 2) (purple), (i = 3,j = 2) (brown), and (i = 3,j = 3) (red) calculated as 
functions of the CO bond distance Ri for i?2 = 2.2 a^ and aoco = 179° along the dashed line 
shown in Fig. |6] 

FIG. 8: Contour maps of the adiabatic state 2^ A": {Ri, R2) plane, with aoco fixed at 179° (a) and 
70° (b), and (i?i,aoco) plane with Rqo fixed at 2.2 oq (c) and 2.6 oq (d). Energy of the dotted 
contour in (a), (c), and (d) is 9.0 eV, in (b) it is 12.5 eV. Contour spacing is 0.25 eV. 

FIG. 9: Contour maps in the (i?i,aoco) plane of the adiabatic states 2^A' (a), 3^yl' (b), l^A" (c), 
and 2>^A" (d). The second CO bond is fixed at 2.2 cq. Energy of the dotted contour is 7.6 eV, and 
the contour spacing is 0.20 eV. 

FIG. 10: Contour maps in the {Ri, R2) plane of the adiabatic states 2^A' (a), 3^A' (b), l^A" (c), 
and 3^A" (d). The angle aoco is fixed at 160°. Energy of the dotted contour is 8.0 eV, and the 
contour spacing is 0.25 eV. 

FIG. 11: CASSCF energies of the electronic states 1 — 5^A' (a,c,e) and 1 — 5^A" (b,d,f) along one 
CO bond distance Rco- The second CO bond is fixed at 2.4 ag: the fixed bond angle is indicated 
in each panel. In all panels, dots and solid lines denote adiabatic energies. In panels (a) and 
(b), lines are color coded according to the ab initio (L^) values as explained in text. Following 
a particular color, one follows a diabatic state. Diabatic assignments using spectroscopic symbols 



are given. Arrows in panel (b) indicate the A/II CIs discussed in Sect. IIIC For A" states, the 
ground electronic state X^A' , shown with gray dotted line, sets the vertical energy scale. 

FIG. 12: Adiabatic CASSCF energies of states 1—5^' (a) and 1—5^A" (b) along the bond angle. 
Fixed CO bond distances are 2.4 oq and 2.5 cq. The two uppermost states in each symmetry block 
are drawn gray. Arrow indicates direction towards the carbene-like OCO minimum in the given 
state. Correlation with the diabatic states at linearity is marked to the right of each panel. 

FIG. 13: Contour maps in the (-Ri,aoco) plane of the adiabatic states 2^A' (a), 3^A' (b), l^A" (c), 
and 3^A" (d). The second CO bond is fixed at 2.6 oq. Energy of the dotted contour is 9.0 eV, and 
the contour spacing is 0.20 eV. 
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FIG. 14: The mixing angle G for the states 2, 3M calculated as a function of one CO bond distance. 



Solid lines and open circles denote G from the regularized diabatic states model of Eqs. (12) — (14). 
Solid squares are the values obtained using quasi-diabatization procedure in MOLPRO. The second 
CO bond is fixed at Ri = 2.2 ao; the bond angle aoco is indicated in each panel. The ab initio 
mixing angle, obtained with the cc-pVQZ basis set, is lifted by 25°. 



FIG. 15: Diabatic potentials (diagonal matrix elements) for the S (green), 11 (blue), and A (red) 
states and the off-diagonal couplings (black solid line) as functions of one CO bond distance. A' 
{A") states are shown in the left (right) panels. The second CO bond is fixed at 2.2 ao; the fixed 
bond angle is indicated in each panel. Black dashed line indicates the ground electronic state. 



FIG. 16: The diabatic potentials Vn' (a), Va' (c), Vw (b), and Va" (d), and the off-diagonal 
coupling elements Fn'A' (e) and Vn"A" (f) in the (i?i, R2) plane with the bond angle fixed at 175°. 
Energy of the dotted contour in (a), (b), (c), and (d) is 7.5 eV, and the contour spacing is 0.25 eV. 
In (e) and (f), the respective numbers are -0.3 eV and 0.05 eV. 



FIG. 17: The diabatic potentials Vn' (a), Va' (c), Vu" (b), and Va" (d), and the off-diagonal 
coupling elements Vq'A' (e) and Vwa" (f) in the (i?i,aoco) plane with the CO bond fixed at 
2.2 cq. Energy of the dotted contour in (a), (b), (c), and (d) is 8.6 eV, and the contour spacing is 
0.20eV. In (e) and (f), the respective numbers are -1.20eV and 0.20eV. 
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